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Abstract 

Members of the TOPLESS gene family emerged recently as key players in gene repression in several mechanisms, 
especially in auxin perception. The TOPLESS genes constitute, in 'higher-plant' genomes, a small multigenic family 
comprising four to 11 members. In this study, this family was investigated in tomato, a model plant for Solanaceae 
species and fleshy fruits. Six open reading frames predicted to encode topless-like proteins (SITPLs) containing the 
canonical domains (LisH, CTLH, and two WD40 repeats) were identified in the tomato genome. Nuclear localization 
was confirmed for all members of the SITPL family with the exception SITPL6, which localized at the cytoplasm and 
was excluded from the nucleus. SITPL genes displayed distinctive expression patterns in different tomato organs, 
with SITPL1 showing the highest levels of transcript accumulation in all tissues tested except in ripening fruit where 
SITPL3 and SITPL4 were the most prominently expressed. To gain insight into the specificity of the different TOPLESS 
paralogues, a protein-protein interaction map between TOPLESS and auxin/indole-3-acetic acid (Aux/IAA) proteins 
was built using a yeast two-hybrid approach. The PPI map enabled the distinction of two patterns: TOPLESS isoforms 
interacting with the majority of Aux/IAA, and isoforms with limited capacity for interaction with these protein partners. 
Interestingly, evolutionary analyses of the TOPLESS gene family revealed that the highly expressed isoforms (SITPL1, 
SITPL3, and SITPL4) corresponded to the three TPL-related genes undergoing the strongest purifying selection, while 
the selection was much weaker for SITPL6, which was expressed at a low level and encoded a protein lacking the 
capacity to interact with Aux/IAAs. 

Keywords: Aux/IAA, auxin signalling, co-repressor, multigenic family, protein-protein interactions. Solarium lycopersicum , 
tomato, TOPLESS. 



Abbreviations: AD, activating domain; ARF, auxin response factor; Aux/IAA, auxin/indoie-S-acetio acid; BD, binding domain; EST, expressed sequence tag; NLS, 
nuclear localization signal; PPI, protein-protein interaction; qRT-PCR, quantitative reverse transcription-PCR; Y2H, yeast two-hybrid; YFR yellow fluorescent protein. 
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Introduction 

It is now well accepted that transcriptional co-repressors play 
crucial roles in a broad range of plant developmental pro- 
cesses (Liu and Karmarkar, 2008; Krogan and Long, 2009). 
In land plants, the Groucho (Gro)/Tupl family of co-repres- 
sors includes TOPLESS/TOPLESS-RELATED (TPL/TPR) 
and LEUNIG/LEUNIG HOMOLOG (LUG/LUH) (Conner 
and Liu, 2000; Kieffer et al, 2006; Long et al, 2006). TPL 
proteins have been shown to be involved in multiple signal- 
ling pathways in higher plants, including hormone-signalling 
pathways (auxin, jasmonic acid, abscisic acid, and ethylene), 
meristem maintenance, floral induction, biotic stress, and 
circadian oscillator mechanism (Liu and Karmarkar, 2008; 
Szemenyei et al, 2008; Pauwels et al, 2010; Zhu et a!., 2010, 
Causier et al., 2012a, b; Wang et al., 2013). 

The first TPL gene was identified in Arabidopsis as respon- 
sible for the semi-dominant tpl-1 embryo development 
mutation resulting in altered polarity, ranging from fused cot- 
yledons to complete replacement of the shoot with a second 
root (Long et al, 2002, 2006). Subsequently, the TPL family 
in Arabidopsis thaliana was found to comprise five members 
that seem to act redundantly (rPL, TPRl, TRP2, TRP3, and 
TRP4). Indeed, a quintuple loss of function, in which all five 
TPL/TPR genes were inactivated by mutation or RNA inter- 
ference, is required to phenocopy the tpl-1 phenotype (Long 
et fl/.,2006). 

It was established that, although TPL proteins are lacking 
a DNA-binding activity, they are incorporated into transcrip- 
tion complexes by interacting with transcription factors to 
repress gene expression in various processes. This inhibition 
of the expression of target genes is mediated by the recruit- 
ment of histone deacetylases into transcription complexes, 
and by changing the chromatin state from active to inactive 
(Long et al, 2006; Liu and Karmarkar, 2008; Krogan and 
Long, 2009; Krogan et al, 2012). Interaction between the 
TPL/TPR co-repressors and transcription factors depends 
on the Lissencephaly (LisH) and the C-terminal to LisH 
Homology (CTLH) domain of TPL (Szemenyei et al, 2008; 
Gallavotti et al, 2010), and on a small conserved protein 
motif found in transcription factors. This motif is known as 
the ethylene response factor-associated amphiphilic repres- 
sion (EAR) domain (Ohta et al, 2001), with the consen- 
sus sequence (L/F)DLN(L/F)xP (Ohta et al, 2001; Hiratsu 
et al, 2004). Recently, the Arabidopsis TPL/TPR interactome 
framework revealed that the TPL co-repressors are able to 
interact with various transcription factors harbouring differ- 
ent repression domains (Causier et al, 2012fl). Among these 
TPL interactants, the transcriptional repressors involved 
in auxin signalling [i.e. auxin/indole-3-acetic acid (Aux/ 
lAA) and auxin response factor (ARE) families] have been 
well documented. In Arabidopsis, the discovery that TPL is 
recruited by Aux/IAA proteins to suppress the expression 
of auxin-responsive genes in the absence of auxin revealed 
a crucial role for TPL in mediating the inhibitory effect of 
Aux/IAA on ARF-regulated transcription (Szemenyei et al , 
2008). Large interactome studies in Arabidopsis identified 20 
of the 29 AtlAA proteins as interacting partners of the TPL/ 



TPRs (Arabidopsis Interactome Mapping Consortium, 2011; 
Causier et al, 2012a). In addition, a large-scale analysis of 
the interaction between Aux/IAA and ARE in the Arabidopsis 
shoot apex revealed that the vast majority of the Aux/IAAs 
interacted with all the ARE activators and showed very lim- 
ited interactions with ARE repressors (Vernoux et al , 201 1). 
However, a recent study showed that repressive ARE pro- 
teins, such as ARE2 and ARE9, can interact directly with 
TPL/TPR proteins, suggesting a mechanism for repression 
implicating TPL/TPR co-repressors in both forms of ARE- 
mediated repression (Causier et al, 2012fl). 

The release in recent years of several plant genome 
sequences has offered the possibility to investigate a large set 
of multigenic families at the genome scale. In this context, 
the tomato genome is of special interest, as (1) tomato has 
emerged as a model plant, for fleshy fruit development, and 
(2) tomato is a reference species for the Solanaceae family and 
also for the taxum of Asterids, particularly as the majority 
of sequenced dicot genomes belongs to Rosids (Sato et al 
2012). It is noteworthy that the structure of several multigenic 
families involved in auxin perception and responses have been 
examined in tomato (Kumar et al, 2011; Ren et al, 2011; 
Wu et al, 2011; Audran-Delalande et al, 2012; Kumar et al, 
2012;Pattison & Catala, 2012; Wu et al, 2012a, b), thus shap- 
ing an exhaustive picture of auxin signalization complemen- 
tary to the Arabidopsis model plant. However, compared with 
the plant model Arabidopsis, the TPL gene family has so far 
been poorly described. 

To characterize fully the molecular biology and evolu- 
tion of the tomato TPL family and to understand its pos- 
sible functions, we identified and characterized six SITPL 
genes. Our analyses focused on the identification, evolution- 
ary relationships, and expression patterns of each member of 
the tomato TPL family. Moreover, we used yeast two-hybrid 
(Y2H) approaches to establish the framework of TPL/ Aux/ 
lAA protein-protein interactions (PPIs). These results will 
provide a framework for further studies to better understand 
the potential functions of TPL proteins in tomato plants, 
especially during the flower and fruit development. 

Materials and methods 

Isolation and cloning of SITPL genes 

The full-length coding sequences of six SlTPLs were amplified 
from mature green fruit cDNA. The primers used were as fol- 
lows: TPLl_attbl: 5'-ATGTCATCTCTCAGTAGAG AGCTT-3' 
and TPLl_attb2: 5'-TCATCTTGGTGCTTGATCGGAGC-3'; 
TPL2_attbl: 5'-ATGTCTTCCTTGAGTAGGGAACTG -3' and 
TPL2_attb2: 5'-TCACCTTGAAGGTGTTTCTGATG-3'; TPL3_ 
attbl: 5'-ATGTCTTCTCTTAGCAGAGAATTG-3' and TPL3_ 
attb2: 5'-TCATCTTTGAACTTGGTCAGCAG-3'; TPL4_attbl: 
5'-ATGACTTCTTTAAGCAGAGAGCTG-3' and TPL4_attb2: 
5'-CTACCTTGATGCTTGATCAAGACC-3';TPL5_attbl:5'-ATG 
AGGCATTTTGATGAAATGGT-3' and TPL5_attb2: 5'-CT 
ACCTTTGAGGTTGATCT GAAT-3'; and TPL6_attbl: 5'-ATGT 
CTCTTAGTAAGGACCTTAT-3' and TPL6_attb2: 5'-CTATATTG 
GTTGCTCAT TGGTAA-3'. 

After amplification, the SITPL genes were cloned into the pDO- 
NOR207 vector using the Gateway method (Invitrogen) and were 
fully sequenced. 
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Subcellular localization of SITPL proteins 

For localization of the SITPL proteins, the SITPL coding sequences 
were cloned using Gateway technology as a C-terminal fusion in 
frame with yellow fluorescent protein (YFP) into the pEarlyGatel04 
vector and expressed under the control of the 35S cauliflower mosaic 
virus promoter. The empty vector pEarleyGatel04 was used as a 
control. Protoplasts were obtained from tobacco (Nicotiana taba- 
cum) suspension-cultured BY-2-cells and transfected according to a 
method described previously (Leclercq et al, 2005). YFP localiza- 
tion by confocal microscopy was performed as described previously 
(Audran-Delalande et al, 2012). 

Expression analysis of SITPL genes 

Total RNA extraction, removal of DNA contamination, cDNA 
generation of eight tomato tissues (root, stem, leaves, bud, flower, 
mature green fruit, breaker fruit, and red fruit), and quantitative 
reverse transcription-PCR (qRT-PCR) were performed accord- 
ing to methods previously described (Audran-Delalande et al., 
2012; Pirrello et al, 2006). The primer sequences were as fol- 
lows: TPLIF: 5'-TGTTCGT TCTAGGAGACTAACCAG-3' and 
5'-TPLlR: AAGACAAACCTTCCCTTC CGA-3'; TPL2F: 5'-CC 
TGTAAATACGCCT CTTGCT-3' and TPL2R: 5'-ACTGGTTGG 
AATGGACTGTG-3'; TPL3F: 5'-CACTTTCTGCTCCAATAA 
CCT-3' and TPL3R: 5'-TCCA TCTGTCAACCCAACTG-3'; TPL4 
F: 5'-CCTTCTAACC CAAGCTCCAG-3' and TPL4R: 5'-AT 
AAACTCCGCCATCAGTA AGTC-3'; TPL5F: 5'-CGTCTATT 
GTAACCCATCCACTC-3'andTPL5R:5'-AGAAGTTACACCAT 
GAGGACCC-3'; and TPL6F: 5'-ACTG GACTAGCATTCTCT 
AACAC-3' and TPL6R: 5'-TTGAATT CCACA CCACTATCTG 
AG-3'. Actin was used as an internal reference. The relative fold 
differences (with SITPL6 as a reference gene) for each sample were 
calculated using the formula 2"'^'^'^'. Three independent RNA isola- 
tions were used for cDNA synthesis and, each cDNA sample was 
subjected to real-time PCR analysis in triplicate. 

Bloinformatic analyses 

SITPL genes were searched using BLAST queries on the Genomic 
(Chromosome v.2.40) and transcript database (cDNA itag 2.4) avail- 
able on the SGN website (http://solgenomics.net/tools/blast/index. 
pi). Exons and introns were deduced from the ITAG 2.3 annotation. 
For S1TPL5 (Solyc07g008040), the 'predicted annotation' missing 
the N-terminal extremity was completed with an additional exon 
(from position 2754093 to 2754173 on SL2.40ch07 chromosome 
7 annotation). Protein domains were first predicted on the pros- 
ite database protein (http://prosite.expasy.org/). Prediction of the 
WD40 segments was refined using the PF00400.27 Pfam Hidden 
Markov Model with an i-value threshold at 0.1. For i-values > 0.1, 
the prediction of WD40 position was deduced from the sequence 
alignment of the different TPL isoforms. 

Nuclear localization signal (NLS) analysis prediction was per- 
formed with 'cNLS Mapper' (http://nls-mapper.iab.keio.ac.jp/cgi- 
bin/NLS_Mapper_form.cgi) (Kosugi et al, 2009). NLS prediction 
scores >5.0 were considered positive. 

Evolutionary analyses 

Phylogenetic analyses and distance matrices were built using the 
MEGA5 package (Tamura et al, 2011). FuO-length amino acid 
sequences were aligned using the ClustalW algorithm. For the 
overall phylogeny, an initial tree encompassing sequences from 
Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, 
Zea mays. Sorghum bicolor, Arabidopsis thaliana. Solarium lyco- 
persicon, Nicotiana benthamiana, Populus trichocarpa, Glycine max 
and Mimulus guttatus was performed using the neighbour-joining 
method. The percentage of replicate trees in which the associ- 
ated taxa clustered together was calculated in the bootstrap test 
(500 replicates). The topology was further confirmed using the 



maximum-likelihood method. Ultimately, a simplified tree was 
performed by limiting the number of genomic sets as the topology 
remained unchanged. Trees were drawn to scale, with branch lengths 
in the same units as those of the evolutionary distances used to infer 
the phylogenetic tree. 

The following genome annotations were used for phyloge- 
netic analyses: Physcomitrella patens (Phypal_l.FilteredModels; 
Reusing et al, 2008); Selaginella moellendorffii (Lycophyte Selmol_ 
GeneModels_FilteredModels3; Banks et al, 2011); A. thaliana 
(TAIRIO; Swarbreck et al, 2008); Populus trichocarpa (Eudicot 
Populus.trichocarpa.v2.0; Tuskan et al, 2006); V. vinifera (12X 
March 2010 release. Glycine max Glymal_pacld; Schmutz et al, 
2010); O. sativa [MSU Rice Genome Annotation (Osal) Release 
6.1; Ouyang et al, 2007]; Z. mays (ZmB73_4a.53_working_transla- 
tions; Schnable et al, 2009); Sorghum bicolor (Sorbil_GeneModels_ 
Sbil_4_aa; Paterson et al, 2009); Solanum lycopersicon 
(ITAG2.3_release; Sato et al, 2012); Brassica rapa (Chiifu-401^2; 
Wang et al, 2011); Eucalyptus grandis (Egrandis_201; http://www. 
jgi.doe.gov/); M. guttatus (Mguttatus_140; http://www.jgi.doe.gov/); 
N. benthamiana (Niben. genome. vO.4.4; Pallas et al, 2012); Solanum 
tuberosum (PGSC_DM_v3.4; Xu et al , 201 1). 

Protein-protein interaction (PPI) assay ofSITPLs and SIIAAs by 
Y2H assay 

Tomato TPL genes were amplified and cloned into the pDBD 
(BD-TPLs) vector (Clontech). Similarly, SlIAA target genes [lAAl 
( JN37943 \),IAA3 ( JN379433), IAA4 ( JN379434), lAA 7 ( JN379435), 
IAA8 (JN379436), IAA9 (JN379437), lAAll (JN379438), IAA12 
(JN379439), IAA14 (JN379441), IAA15 (JN379442), IAA16 
(JN379443), IAA17 (JN379444), IAA19 (JN379445), IAA22 
(JN379447), IAA26 (JN379449), IAA27 (JN379450) and IAA29 
(JN379451)] were inserted in pGAD (AD-IAAs) vectors (Clontech). 
Diploids were selected on medium lacking Trp and Leu, and interac- 
tions were validated by the use of HISS and ADE2 reporter genes 
on medium lacking Trp, Leu, His, and Ade. Manipulation and 
analysis of the Y2H assay followed the manufacturer's instructions 
(Clontech Yeast Protocols Handbook), and all experiments were 
repeated three times independently. For SlTPLl genes lacking LisH, 
the coding sequence was truncated at nucleotide position +112. 



Results 

Identification and cloning ofJPL-related genes in the 
tomato genome 

An in silico search was performed on the tomato genome 
and transcript databases (http://www.solgenomics.net/) using 
Arabidopsis TPL and TPR sequences as queries for BLAST 
searches. While the initial screen identified nine ORFs pre- 
dicted to encode putative TPL-like proteins (SlTPLs), only 
six corresponded to full-length proteins containing all canon- 
ical motifs that define the TPL proteins (Table 1). The full- 
length cDNA of the six SlTPLs was further confirmed by 
RT-PCR amplification, indicating that the corresponding 
coding sequences range from 3396 to 3669 bp with deduced 
protein sizes ranging from 1131 to 1222 aa (Table 1). 

Structural analysis of the six SITPL genes showed that they 
displayed similar numbers of introns (23-25) and exons (24— 
26), except for SITPL6, which was longer than the other TPL 
members (Table 1). Pairwise comparison of the six SITPL pro- 
tein sequences showed that the percentage identity among fam- 
ily members ranged from 44 to 75%. Protein domain searches 
in the Pfam database (http://pfam.sanger.ac.uk/) indicated 
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that all SlTPLs displayed the conserved LisH and CTLH 
domains and had two domains containing several WD40 
repeats: WD40-repeat-l and WD40-repeat-2 with seven and 
five WD40 segments, respectively (Fig. 1 and Supplementary 
Fig. SI available at JXB online). The CTLH domain and the 
WD40-repeat-l were separated by a proline-rich region. 

The tomato TPLs were distributed on four chromosomes: 
two SlTPLs (Solyc03gll6750 and Solyc3gll7360) on chro- 
mosome 3, two (Solyc08g076030.2.1 and Solyc08g029050.2.1) 
on chromosome 8, one (Solyc01gl00050.2.1) on chromosome 
1 and one (Solyc07g008040.2.1) on chromosome 7. There 
were three additional truncated TPL sequences lacking the 
LisH and CTLH domains, with two located on chromosome 
3 (Solyc03gll7370 and Solyc03gll7410) and one on chromo- 
some 1 (Solyc05g016070). 

The number of 'full-length' TPL genes in tomato fell in 
the range found in other plant genomes, which varies in 
angiosperms from four members in monocots to 11 members 
in soybean (Fig. 2). It is noteworthy that a high number of 
isoforms is often observed in organisms having undergone 
recent whole-genome duplication or polyploidization events 
(e.g. G. max, N. benthamiana and B. rapa). 



SITPL nomenclature and phylogenetic analyses 

To adopt a nomenclature consensual with that of Arabidopsis 
TPL and TPR proteins, we carried out phylogenetic analy- 
ses on different TPL-like proteins or cDNAs from different 
plant sequenced genomes comprising moss, fern, and vari- 
ous angiosperm sequences (see Materials and methods). The 
phylogenetic trees (Fig. 3 A) allowed the individualization of 
four branches. Three branches looked well defined in all dicot 
plants: the first branch contained AtTPL, AtTPLl, AtTPR4, 
Solyc3gll7360.2.1 (named SlTPLl), Solyc03gll7360.2.1 
(named S1TPL4), and Solyc07g008040.2.1 (named S1TPL5); 
the second branch, absent in Arabidopsis yet present in 
Eucalyptus (Eucgr.K00093.1|PACid:23601479) and grapes 
(GSVIVT01024440001), contained Solyc08g076030 (named 
S1TPL2), rice ASPl protein, and moss or lycophyte TPL- 
like proteins; and the third branch contained AtTPR2, 
AtTPR3, and Solyc01gl00050.2.1 (named S1TPL3). Lastly, 
Solyc08g029050.2.1 (named S1TPL6) appeared as an out- 
group branch in the phylogenetic tree (Fig. 3A). The robust- 
ness of the tree topology was assessed either with a bootstrap 
test (Fig. 3A) or by changing the number of genomes used 



Table 1 . Main structural features of the tomato SITPL family members 



Nomenclature 


Gene 






Predicted protein 


Domains 






SlTPLs 


iTAG Gene ID 


Exons 


Introns 


Length (aa) 


MW (kDa) 


LisH 


CTLH 


WD-40 repeats 


SITPL1 


Solyc03g1 17360.2.1 


25 


24 


1131 


124.676 


4-36 


34-92 


411-632/832/957 


SITPL2 


Solyc08g076030.2.1 


25 


24 


1136 


124.60 


4-36 


34-92 


341-668/834-959 


SITPL3 


SolycOlgl 00050.2.1 


25 


24 


1132 


124.676 


4-36 


34-92 


343-669/871-955 


SITPL4 


Solyc03g1 16750.2.1 


26 


25 


1133 


124.318 


4-36 


34-92 


413-634/839-964 


SITPL5 


Solyc07g008040.2.1 


24 


23 


1134 


124.82 


4-36 


34-92 


398-639/881/965 


SITPL6 


Solyc08g029050.2.1 


33 


32 


1222 


134.181 


3-35 


33-91 


531-664/934-1060 



iinnnoflEnniiflnnflnfflflnnfflfliflflfa 

SI-tPLY ~-' '•' ~-' '•' ~-' '■' ■•' '•' '•' '•' '•' ■•' ■•' ■■' '•' '•' 

SI-TPL3 ■•' -•- -■- -•- '•■ '•• ■•' '■• -■' ■•• ■•' '•• -.' -■' ■•' -■' -•- 

SI-TPL4' ■■' ■•' ■•' ~-' '•' ~-' ■•' ~-' ■■' ■•' 

SP=TPL5 "•' ■■' '■' '■' '•' ■■' '■' ■■' ~-' ■■' ■•' ■■• ■■' 

Sl^tPlJ6 '■' ■'' ■'' '■' ■-' '■' '■■ ■'' '-' '■' ■■' ''' ''' ■'' 



exon 


■■„ intron „ 


LISH 




WD40-repeat1 



Fig. 1. Gene structure of the six tomato TPL genes. Grey boxes represent exons, dotted lines represent introns, tine red box is tlie LisH 
domain, tlie blue box is the CTLH domain, the light green boxes are the WD40-repeat 1 and the dark green boxes are the WD40-repeat 
2. The figure was produced using FancyGene software (http://bio.ieo.eu/fancygene/). (This figure is available in colour at JXB online.) 
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Fig. 2. Inventory of TPL genes in different plant genomes. Only TPL genes containing the four canonical domains (LisH, CTLH and two 
WD40 repreats) were considered. The major taxons are shown below. 



in the phylogeny and the portion of the ahgned sequence 
(N-terminal, C-terminal, or conserved domains) or the clus- 
tering method (neighbour-joining or maximum-HkeHhood 
method). The vast majority of the nodes presented in Fig. 3A 
remained unchanged. 

To understand further the TPL phylogeny, and notably 
to characterize the S1TPL6 outgroup, the presence of TPL 
'orthologues' was investigated in Asterid genomes belong- 
ing either to the Solanaceae family {Solarium tuberosum and 
N. benthamiana) or to the Lamiales order (M guttatus). An 
S1TPL6 homologue was found in all Asterids, supporting the 
view that S1TPL6 homologues form a distinct clade (Fig. 3B). 
Within this S1TPL6 clade, the length of the branches sug- 
gested that these isoforms had evolved faster than other 
TPLs. This observation was supported by sequence diver- 
gences: the amino acid substitution rates calculated within 
the Solanaceae orthology groups varied from 2.6 to 6.3% 
for SlTPLl-5 and reaching 22.7% for the S1TPL6 (Table 2). 
Moreover, a neutrality test (dS/dN values) calculated on 
Solanaceae orthologues suggested that the purifying selection 
exerted by evolution on the SITPL6 family is much weaker 
than the selection pressure exerted on other TPL genes. 

Subcellular localization ofSITPLs 

The subcellular localization of the SITPL proteins was 
assessed by a transient expression assay in tobacco pro- 
toplasts using a translational fusion between each of the 
SITPL proteins and YFP. Microscopy analysis showed that 



S1TPL1-5-YFP fusion proteins locaHzed exclusively to the 
nucleus (Fig. 4) whereas S1TPL6 was localized at the cyto- 
plasm and excluded from the nucleus. This result is in agree- 
ment with the in silica prediction of a conserved NLS for the 
five nuclear SlTPLl-5 proteins, while S1TPL6 NLS scores 
were below the 5.0 threshold value (Supplementary Table SI 
available at JXB online). Altogether, the nuclear localization 
of the majority of SlTPLs was consistent with their putative 
role in transcriptional regulation activity. 

Expression analyses 

In order to study the spatio-temporal expression pattern of the 
six SITPL genes, qRT-PCR was performed on eight different 
plant tissues and organs. Three SITPL members (SlTPLl, 
SITPL3, and SITPL4) displayed significantly higher levels 
of expression than the three remaining paralogues. SlTPLl 
and SITPL4 were found to be highly expressed in flowers and 
vegetative tissues (roots, stems, and leaves) and in developing 
flowers (buds and during anthesis) but with reduced expres- 
sion in ripening fruit, while SITPL3 expression remained con- 
stant and high during fruit ripening (Fig. 5). This preferential 
expression of SlTPLl, SITPL3, and SITPL4 is coherent with 
their estimated expression in two public databases (RNAseq 
database: http://ted.bti.cornell.edu and ESTs database: http:// 
solgenomics.net/). Although less expressed, SlTPLl was 
found preferentially in leaves and developing flowers; the 
levels of SITPL5 transcripts were low in all tissues; SITPL6 
expression was restricted to roots and stems (Fig. 5b). 
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Fig. 3. Phylogenetic trees of some plant and tomato TPL proteins. (A) Representative pliylogenetic tree of TPL proteins from land 
plants: moss {P. patens, PHYPADRAFT_xxx), lycophyte {Selaginella moellendorffii, SELI\/IODRAFT_xxx), rice (LOC_Os-xxx), tomato 
(red boxes) and Arabidopsis (green boxes). The coloured brackets emphasize the main branches conserved among angiosperms. 
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among Asterid and Solanaceaous species. The tree was built using sequences from four genomes: Solanum lycopersicon, Solanum 
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Table 2. Evolutionary features of TOPLESS-related genes in Soianaceous species 

Mean distance was expressed as the proportion of amino acid or nucleic acids positions different after pairwise alignment. dS/dN values were 
calculated using the codon-based test of purifying selection performed on each pair of orthologous sequences from Solanum lycopersicon and 
Solanum tuberosum. The variance of the difference was computed using the bootstrap method (500 replicates). Analyses were conducted using 
the Nei and Gojobori (1986) method. 



SITPL1 SITPL2 SITPL3 SITPL4 SITPL5 SITPL6 

Mean distance {Solanum/Nicotiana) Amino acids 0.026 0.041 0.032 0.029 0.063 0.227 

Nucleic acids 0.055 0.050 0.054 0.054 0.067 0.154 

Neutrality test (So/amm) dS/dN 7.08 6.66 6.98 7.62 6.19 3.645 



Examination of PPIs in tlie frameworl< of auxin 
mediation 

The differential expression of SITPL genes evokes the criti- 
cal question of functional redundancy within the TPL family. 
In a recent paper, Causier et al. (2012(7) compared the PPI 
patterns of different Arabidopsis TPL proteins using a high- 
throughput Y2H screen both on a whole-plant and on a tran- 
scription factor library. In the present work, we focused on 
the interactions with the Aux/IAA family by performing an 
exhaustive targeted analysis of Aux/IAA-TPL interactions. 
The six SITPL proteins were fused to a binding domain (BD) 



and used as bait in a Y2H test with 17 different SlIAA pro- 
teins fused to an activating domain (AD). After monitoring 
the yeast growth on two auxotroph selective media, two pat- 
terns of TPL could clearly be defined (Fig. 6A, B): SlTPLl, 
S1TPL2, S1TPL4, and S1TPL5 interacted with the majority 
of SlIAAs and grew in all the selective media, and S1TPL3 
and S1TPL6 exhibited only limited growth when co-expressed 
with Aux/IAA-AD fusion proteins. Contrary to other 
SlIAAs, S1IAA29 failed to show interaction with any of the 
SlTPLs. With the exception of S1IAA12 and S1IAA15, the 
Aux/IAAs did not harbour any obvious specificity towards 
the 'TPL' clade (SlTPLl, S1TPL4, and S1TPL5), sharing high 



Tomato Topless family | 1019 



similarity with AtTPL. In addition, S1TPL2, which belongs 
to a distinct clade of SlTPLs (1,4, and 5), also exhibited a 
broad capacity to interact with the majority of SlIAAs. As 
a control, we performed a Y2H test with truncated SlTPLl 
or S1TPL5 (ALisH-TPL) (Fig. 6C) lacking the LisH domain 
shown previously to be essential for TPL-WUS or TPL-Aux/ 
I AA interactions (Kieffer et al. , 2006; Szemenyei et al , 2008). 
Contrary to all SlTPLs BD fusions assayed, a complete lack 
of growth was observed when co-expressing BD-ALisH-TPL 
proteins with BD-Aux/IAAs (Fig. 6B). 



Discussion 

The present study addressed the structural, evolutionary, and 
functional features of the tomato TPL family. TPL proteins 
have been primarily defined as a major component of the 
auxin transduction and response pathway, but the present 
data sustain the hypothesis of a functional diversification of 
these regulatory proteins. While mainly focusing on the TPL 
family in tomato, a plant model for Solanaceae and fleshy fruit 
research, the data also addressed the comparative features of 
this gene family within Plant kingdom at the evolutionary 
level, shedding new light on their functional diversification. 

The structure of the SITPL family is representative of 
that found in angiosperms where these proteins belong to a 
small multigenic family comprising five to 11 members. In 
the tomato, six full-length SITPL genes were identified, as 
well as additional three pseudogenes with incomplete coding 
sequences. Among the six SITPL genes, five were highly con- 
served {SlTPLl-5), while the last gene {SITPL6) was more 
distant. With the exception of poplar genomes and genomes 
having undergone recent polyploidization (i.e. soybean, 
B. rapa, and N. benthamiana), the number of TPL isoforms 
ranges from four to six members, suggesting that the number 
of genes remains stable in this family and that usually, after 
a whole-genome duplication event, duplicated copies of TPL 



genes are not retained. The phylogenetic analysis of TPL 
genes enabled the distinction of three major clades gather- 
ing homologues in the majority of angiosperm genomes. The 
last clade, containing the distant SITPL6, displays only clear 
homologues in closely related taxa (Asterids). Interestingly, 
highly diverging sequences of TPL-related proteins have also 
been found in other genomes such as the AtTPR-Uke gene 
(At2g25420; Causier et al, 2012b) and in poplar, but no clear 
relationship could be established with SITPL6. Contrary to 
angiosperm TPL proteins, TPL from Physcomitrella patens 
and Selaginella moellendorffii clustered in a same branch, 
indicating the existence of ancestral divergences occurring 
before angiosperm radiation. 

The functionality of SITPL genes was addressed through 
three approaches: expression analysis, subcellular locali- 
zation, and establishment of an interaction map between 
SITPL and SlIAA proteins. The expression patterns of differ- 
ent SlTPLs revealed the tissue specificity of various isoforms 
and suggested a functional specialization of SITPL isoforms. 
For example, SlTPLl is highly expressed in vegetative organs 
(stems, roots) and flowers, while the expression of SITPL3 
and SITPL4 prevails in fruit. Moreover, the overall inten- 
sity of gene expression evaluated by qPCR demonstrated 
a distinction between a group of three isoforms {SlTPLl, 
SITPL3, and SITPL4) that are highly expressed, SITPL2, 
which is moderately expressed in the leaves and flowers, and 
third group made of two isoforms (SITPL5 and SITPL6) that 
displayed very low levels of expression. In agreement with our 
data, the prevalence of SlTPLl, SITPL3, and SITPL4 tran- 
scripts was also observed in expressed sequence tag (EST) and 
RNAseq expression databases (http://ted.bti.cornell.edu), 
whereas the expression of SITPL6 was again found to be very 
low (no EST and few RNAseq reads). Interestingly, the over- 
all expression level negatively correlated with the amino acid 
substitution rate. Indeed, after defining orthology groups 
among Solanaceous TPLs, we found that the highly expressed 
isoforms {SlTPLl, SITPL3, and SITPL4J showed the highest 




Fig. 4. Subcellular localization of tomato TPL proteins. SITPL-YFP fusion proteins were transiently expressed in BY-2 tobacco 
protoplasts and subcellular localization was analysed by confocal laser-scanning microscopy The merged pictures of the yellow 
fluorescence channel (left panels) and the corresponding bright field (middle panels) are shown (right panels). The empty vector 
pEarleyGatel 04 was used as a control. Bar, 1 0 )im. (This figure is available in colour at JXB online.) 



1020 I Hao etal. 



B 



— 150 
Z100 

E 50 
S 0 



Root 



5 <?^<?^<?^<?'^<?'^^<?^ 

l£ <^ <^ 



fc40 
I 30 
K 20 
S 10 
S 0 
IS -10 



stem 



^ <^ 



— 600 
I 400 
E 200 

I " 

1-200 



Leaves 






2150 

fioo -\ 



I 50 



SITPL1 



■ <^ 1^ ;« 



z 4 

K 

E 2 

I 0 



SITPL2 



■ill.. 



' <? * i> 

a? xs^ „< 



Bud (Flower) 






~ Anthesis (Flower) 



£800 
< 600 
K 400 
E 200 
I 0 



vT- v'b .b' .fe 
cj^ 



Mature Green Fruit 



40 
30 
20 
10 
0 

-10 J\ ^-b .b .6 

^ ^ ^ 



Lk 



SITPL3 



SITPL4 



6 1 
4 - 
2 

0 - 



<;3o 

1io I- I 



Breaker fruit 



< 80 

K 60 

i '° 

■M 20 

■5 0 



.1 v'b vt' v*© 
^ <^ £^ 



Red fruit 



Z 200 T 

1^ 150 ■ 
Li full 



<^ <^ 



a 2 

< 1,5 

1 1 

s 0,5 

« 0 



SITPL5 



120 

<i5 ^ 
1 10 

■I 0 



SITPL6 
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amino acid sequence conservation (<3.2% difference within 
Solanaceous sequences), while sequences were less conserved 
within the S1TPL6 orthology group (22.7% difference within 
Solanaceous sequences). The moderately expressed SITPL2 
and SITPL5 displayed intermediate substitution rates (4 and 
6% differences, respectively). This correlation was also sup- 
ported by a neutrality test ((iS/JA^ values) performed between 
potato and tomato pairs of othologues. The high substitution 
rate within the SITPL6 orthology group was interpreted as an 
indication that the SITPL6 subfamily undergoes a reduced 
purifying selection. By contrast, broadly expressed SITPL 
isoforms are under a stronger purifying selection. Such a 
correlation between gene expression level and amino acid 
substitution rate has already been observed in genome-wide 
comparisons of expression patterns and protein evolution in 
Ambidopsis-related plants and in the Poaceae family (Wright 
et al, 2004; Slotte et al. , 201 1; Davidson et al. , 2012). Indeed, 
this correlation is consistent with A. thaliana expression data 
(AtGenExpress), At-TPL being expressed more than other 
AtTPRs and At TPL orthologues remaining highly conserved 
either in Arabidopsis lyrata or in B. rapa. 

The subcellular localization established a second discrimi- 
nation criterion among SlTPLs. YFP fusion proteins of 
SITPL 1-5 isoforms all migrated exclusively to the nucleus, 
as observed with other TPL proteins from Arabidopsis (Long 
et al, 2006), maize (Gallavotti et al , 2010), and rice (Yoshida 
et al, 2012). By contrast, the S1TPL6-YFP fusion protein 
displayed a divergent subcellular targeting, this isoform being 
targeted to the cytosol. This divergent localization is in line 



with the lower scores calculated by the NLS prediction tool 
for S1TPL6. This observation, in addition to the low expres- 
sion level and the high substitution rate, supports the view of 
either a partial loss of functionality or divergent functionality 
regarding S1TPL6. 

The first established function of TPL proteins is related to 
their role in auxin signalling via interaction with Aux/IAA 
partners (Szemenyei et al, 2008). To check whether this role 
is conserved among all SlTPLs isoforms and gain insight 
on either functional redundancy or potential functional 
diversification among family members in tomato, a compre- 
hensive PPI study was carried out between all SlTPLs and 
SlIAA members using a Y2H screen. This targeted interac- 
tome study revealed two distinct patterns of interaction for 
tomato TPLs: four isoforms (SlTPLl, S1TPL2, S1TPL4, and 
S1TPL5) displayed a broad capacity for interaction with the 
majority of SlIAAs, and the remaining two isoforms (S1TPL3 
and S1TPL6) showing a more restricted interaction capacity. 
It is noteworthy that a large number of SlIAAs showed posi- 
tive interaction with SlTPLs, consistent with the outcome of 
Y2H screens performed in Arabidopsis where 20 out of the 29 
AtAux/IAAs were able to interact with AtTPLs (Szemenyei 
et al, 2008; Arabidopsis Interactome Mapping Consortium, 
2011; Causier et al, 2012fl). Interestingly, neither S1IAA29 
nor its Arabidopsis homologue At IAA29 (AT4G32280.1) 
interacted with TPL proteins, although S1IAA29 exhib- 
its a repressor activity (Audran-Delalande et al, 2012). On 
the other hand, the limited interaction capacity displayed 
by S1TPL6 adds another distinctive feature to this isoform. 
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Fig. 6. PPI maps between SITPLs and SIIAAs established by a Y2H screen. (A) Yeast growth of co-transformed BD-TPLs and AD-IAAs. 
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which has aheady diverged from other family members by 
its low expression level, high amino acid substitution rate 
and different subcellular localization. Altogether, the cumu- 
lative distinctive features support the idea that S1TPL6 has 
partially lost its ancestral function and may have gained new 
functionality 

In previous Y2H screens performed in Arabidopsis by 
Causier et al (2012a), AtTPR3 and AtTPR2, closely related 
to S1TPL3, both displayed the capacity to interact with vari- 
ous Aux/IAA proteins. However, a closer look at the inter- 
action map pubhshed by Causier et al. (2012o) could also 
suggest differences in specificity between AtTPL and AtTPR2 
or AtTPR3, with the two latter notably interacting with part- 
ners displaying partial repression domains. Such hypothesis 
opens the possibihty that At-TPR2, At-TPR3 and the closely 
related S1TPL3 display a specialization alternative to auxin 
signalling. The development of quantitative PPI methods 
such as Forster resonance energy transfer or surface plasmon 



resonance may provide deeper insight on discriminating 
interaction features among various TPL isoforms. 

Functional redundancy among Arabidopsis TPL family 
members is supported by the absence of obvious phenotypes 
in single loss-of-function mutants of AtTPLITPR genes and 
by the requirement for downregulation of all five AtTPL- 
TPRs in order to phenocopy the dominant mutation tpl-1 
(Long et al., 2006). However, this assumption is in contrast 
to the situation prevailing in rice and maize, where genetic 
evidence seems to support a more specialized functionality 
for TPL genes. Thus, in rice (Yoshida et al., 2012), a single 
recessive mutation in Aspl, a TPL-like gene close to SITPL2, 
exhibited several pleiotropic phenotypes, such as altered 
phyllotaxy and spikelet morphology While these phenotypes 
suggest a close association of Aspl with auxin action, they 
clearly reveal that the specialization of TPL-related pro- 
teins in some organisms can differ from that in Arabidopsis. 
Further evidence sustaining a diversified function for TPL 
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proteins is provided by maize rel2 mutants affected in a TPL- 
like gene closely related to SITPL3 and AtTPRS (Gallavotti 
et al, 2010). A better clarification of the putative specialized 
functionality among tomato TPLs might be addressed by a 
reverse genetics approach. Simultaneous downregulation of 
SlTPLl and SITPL4 would uncover the importance of the 
TPL family in vegetative development and auxin action. 
Likewise, specific downregulation of SITPL3 would be of 
particular interest to unravel the role of TPL co-repressors in 
flower and fruit biology. 

Altogether, these data shed new light on structural, evolu- 
tionary, and some functional features of the tomato TPL gene 
family that suggest functional diversification of these regulatory 
proteins. Of particular interest, the setup of a comprehensive 
TPL-Aux/IAA interaction map and the differential subcellu- 
lar targeting of some SlTPLs proteins would provide important 
clues towards designing appropriate strategies for the elucida- 
tion of both redundant and specific roles of TPL genes. 

Supplementary data 

Supplementary data are available at JXB online. 

Supplementary Fig. SI. Multiple sequence alignment of 
full-length SITPL proteins. 

Supplementary Table SI. NLS prediction scores computed 
with cNLS Mapper (Kosugi et ah, 2009). 
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